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Abstract 

We prove short-time existence of smooth solutions for a class of 
nonlinear, and in general spatially nonlocal, Hamiltonian evolution 
equations that describe the self-interaction of weakly nonlinear scale- 
invariant waves. These equations include ones that describe weakly 
nonlinear hyperbolic surface waves, such as nonlinear Rayleigh waves 
in elasticity. 



1 Introduction 

In this paper we prove the short-time existence of solutions of a class of 
quadratically nonlinear evolution equations in one space-dimension. The 
equations are, in general, spatially nonlocal; an example is the following 
singular integro-differential equation for (p(x, t) 

^ t + ^U[ij 2 ] xx + ^ xx = 0, (1.1) 

where H denotes the Hilbert transform with respect to x and ip = H [</?]. 
According to the result we prove here, the initial value problem for this 
equation has a unique short-time solution in the Sobolev space H s when 
s > 5/2. 
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This class of equations arises from a Hamiltonian description |33| of 
weakly nonlinear, translation-invariant, scale-invariant waves 0. By 'scale- 
invariant' we mean that the wave motion is invariant under space-time 
rescalings x i— > rjx, t » fyi, for all 77 > 0, with the dependent field vari- 
ables scaling according to their appropriate space-time dimensions. As a 
result, the wave motion does not depend on any intrinsic length or time 
scales, and any space-time parameters on which it does depend must have 
the dimension of velocity. 1 One consequence of this scale- invariance is that 
the wave motion is non-dispersive. 

A specific physical motivation for studying this class of equations comes 
from quasi-linear hyperbolic surface waves, such as Rayleigh waves in elas- 
ticity. 2 Surface waves that satisfy half-space or planar interface problems 
for systems of conservation laws (without source terms) are scale-invariant 
because neither the PDE nor the geometry define an intrinsic length scale. 

We refer to surface waves, like Rayleigh waves, that decay exponen- 
tially away from the boundary on which they propagate, as 'genuine' surface 
waves, in contrast to 'radiative' or 'leaky' surface waves that propagate bulk 
waves into the interior. Weakly nonlinear 'genuine' surface waves are de- 
scribed asymptotically by nonlocal equations, of which is an example, 
whereas purely 'radiative' surface waves are typically described by standard 
local PDEs (see [3] for an example) . One of our main goals is to prove short- 
time existence for such nonlocal asymptotic surface wave equations, and it 
is useful to put them in the context of a broader Hamiltonian framework. 

These surface wave equations are of interest, in particular, because they 
describe the formation of singularities in quasi-linear hyperbolic waves on 
the boundary of a half-space, a mechanism that differs from the familiar one 
of shock formation in the interior of a domain. Their study should also shed 
light on the nonlinear well-posedness of IBVPs for hyperbolic PDEs, and 
on the stability of shocks and other interfaces, when the uniform Lopatinski 
condition fails as a result of the presence of 'genuine' surface waves. (See pQ, 
[2] for further discussion, and [2Ej for a local existence result for nonlinear 
boundary value problems in elasticity that permit Rayleigh waves.) 

The canonical spectral form of the wave equation we study here is given 
in (|2.4|) . This equation describes the evolution of a wavefield subject to 
three-wave resonant interactions between wavenumbers {k, m, n} such that 

1 Wave motions, such as gravity water waves on deep water, that depend only on an 
acceleration rather than a velocity are not scale-invariant in the sense we use the term 
here. 

2 The term 'surface wave' is also used to refer to water waves, whose properties are 
different from the non- dispersive surface waves we discuss here. 
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k + m+n = 0. These three- wave resonances always occur for scale-invariant 
waves, since they are non-dispersive. Equation (|2,4|) depends on a given ker- 
nel, or interaction coefficient, T : 1? — > C. Provided that T(k, m, n) is not 
identically zero when k + m + n = 0, we expect that the dominant weakly 
nonlinear effects are those given in (|2.4|) ; thus, this equation provides a gen- 
eral description of nonlinear, scale-invariant waves in one space dimension 
when higher-order resonances can be neglected. 

The Hamiltonian formalism used to derive Q2.4|) does not imply that the 
resulting equations are well-posed. Short-time existence is not obvious, and 
suitable bounds on the interaction coefficients appear to be required. (See 
|13j for some examples of related nonlocal, but non-Hamiltonian, evolution 
equations that are ill-posed, and the discussion of compacton-type equations 
in Section 13.41 ) Consequently, the result we prove here provides a useful 
validation of the corresponding equations as consistent models for the self- 
interaction of nonlinear waves. 

Solutions of these equations typically form singularities in finite time, 
so that, in general, a short-time existence result for smooth solutions is the 
best one can hope for. Singularity formation and the global existence of 
appropriate weak solutions are interesting questions, but we do not address 
them here. Unlike local existence, global existence presumably depends on 
the specific kernel, making a general global existence theory unlikely. 

In proving our existence result, we find it convenient to study a non- 
canonical form of (|2.4|) . given in (|2,lll) . which may be regarded as a nonlocal 
generalization of the inviscid Burgers equation. The proof uses a standard 
Galerkin procedure, although we work in Fourier space because of the spa- 
tially nonlocal structure of the evolution equation. 

The Hamiltonian structure of Q2.4[) imposes certain reality and symmetry 
conditions on the kernel T, given in (|2,5|) . In particular, it implies that T 
is symmetric under under cyclic permutations of its arguments, and this 
property plays a crucial role in our short-time existence proof. 

Scale-invariance further implies that T is a homogeneous function whose 
degree v is a parameter that characterizes the space-time scaling invariance 
of the equation. This parameter takes the value 3/2 for bulk waves and 2 
for surface waves, but our analysis applies to kernels of arbitrary degree. 

It is convenient to treat the cases v > 3/2 and v < 3/2 separately. For 
v > 3/2, we assume that the interaction coefficients satisfy the bound in 
(|4.9|) . This condition appears to be necessary to develop a short-time H s - 
existence theory in this case, and it is satisfied by the kernels that arise in 
surface wave problems. Our existence result is stated in Theorem ^ and 
the main estimate is given in Proposition El One quantity of interest is the 
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number of square-integrable spatial derivatives required to prove short-time 
existence. As we show, short-time existence holds in H s for the noncanonical 
equation (|2.11f) when s > v. 

An existence result for v < 3/2 under the assumption that the interaction 
coefficients satisfy the bound in (|4.3U|) — a weaker condition than (|4.9j) - 
is given in Theorem 

In Section we write out the general equation in various forms (canon- 
ical and noncanonical, spatial and spectral). In Section |HJ we give some 
examples that arise from specific choices of the kernel T. These include the 
inviscid Burgers equation, the Hunter-Saxton equation, asymptotic equa- 
tions for nonlinear surface waves, and equations related to compacton equa- 
tions. In Section |IJ we prove the local existence result, and in Section |S] we 
state an inequality used in the proof. 

2 The evolution equation 

We consider an equation for a scalar-valued function p(x, t) of the form 



Here, H denotes the Hilbert transform with respect to the one-dimensional 
spatial variable x, and 6(93, ip) is a quadratically nonlinear function. The 
skew-adjoint Hilbert transform plays the role of a Hamiltonian operator, 
and b(<p, ip) is the variational derivative of a cubically nonlinear Hamiltonian 
functional. 

The specific structure of the equation is most easily described in Fourier 
space. We study spatially periodic solutions for definiteness, so that 



where T is the circle of length 2ir, and / = (— i*) is a time interval. Most 
of the results extend in a straightforward way to equations defined on M 
instead of T. 

We denote the Fourier coefficients oft^by^iZx/— >C, where 



tpt + H [b(p, tp)] = 0. 



(2.2) 



(p : T x I -> R 



00 



k=— 00 




4 



The Hilbert transform is defined by 

H[(p](x,t) = -^p.v. J^cot (^—^J <p{y,t)dy, 
H[<p](k, i) = isgn k <p(k, t), 

where 

( 1 k > 0, 
sgn k = < k = 0, 
[ -1 k < 0. 

The bilinear form 6 is defined in Fourier space by 

oo 

b(Jp^jj)(k,t) = ^2 T(-k,k-n,n)(p(k-n,t)$(n,t), (2.3) 

71= — OO 

where T : Z 3 — > C is a given kernel. Only the values of T(k, m, n) with 
k + m + n = appear in 1)2. 3 jl . so we could omit one of the arguments, but 
it is convenient to include all of them to exhibit the Hamiltonian symmetry 
property stated below. It then follows from (|2.2[) and (|2.3|) that (p{k, t) 
satisfies the evolution equation 

OO 

(pt{k, t) + isgn k T{—k,k — n,n)(p(k — n,t)(p(n,t) = 0. (2.4) 

n=— oo 

We assume that T : Z 3 — > C satisfies the following reality, symmetry, 
Hamiltonian, and scaling conditions: 

T(k, m, n) = T*(—k, —m, —n), 

T(k, m, n) = T(k, n, m), (2-5) 
T(k, m, n) = T(m, n, k), 

T(r]k,7]m,rjn) = r) v T(k,m,n) for all rj > 0. 

Here, a star denotes the complex-conjugate. The symmetry of T(k, m, n) 
under cyclic permutations of {k, m, n} gives ()2.4|) a Hamiltonian structure. 
For simplicity, we also assume that 

T(0, m, n) = T(k, 0, n) = T(k, m, 0) = 0, 

so that there are no mean- field interactions. 
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The homogeneity of T of degree v corresponds to an invar iance of Q2.2|) 
under rescalings 

t i— » rjt, ryx, ip i — > rf^ip. (2-6) 

Dimensional analysis 2 shows that the parameter 1/ is related to the co- 
dimension (n — (f) of the waves by 

(2.7) 

Here, d is the dimension of the wavenumber space, and n is the spatial 
dimension of the parameter p that sets the mass dimension of the wave 
motion; for example, if this parameter is a density, then n is the dimension 
of the space in which the waves propagate. It follows that v = 3/2 for bulk 
waves with d = n, and v = 2 for surface waves with d = n — 1. 

The dimensional analysis also shows that if (p is a canonical variable, 
then ip/p has dimension L^ _1 / 2 T -1 / 2 , where L and T denote space and 
time dimensions, respectively. This explains the scaling of (p by r] u ~ 1 in 



2.1 Hamiltonian structure 

For simplicity, we consider functions whose spatial mean is zero, so that 
<p(0,t) = 0. (This involves no loss of generality, since p(0,t) is constant in 
time.) We define the Hamiltonian 

oo 

7i{p,(p*) = {T(—m — n,m,n)p* (m + n,t)tp(m,t)tp(n,t) 

m,n=l 

+ T*(—m — n, m, n)<p*(m, t)<p*(n, t)<p(m + n, t)} . 

The complex canonical form of Hamilton's equations, 

81~L 

m(k,t) = —(k,t) forfceN, 

then gives 

fc-i 

ipt(k, t) = T(—k, k — n, n)p(k — n, t)tp(n, t) 

n=l 

oo 



+2^T*(-fc - n,k,n)cp(k + n,t)tp*(n,t) 



n=l 
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The first term on the right-hand side of this equation describes the upward 
transfer of energy by nonlinear interactions to higher wavenumbers, while 
the second term describes the reverse downward transfer. These equations 
may be rewritten as (|2.4|) for k € 7L by use of the properties of T and 
the fact that ip*(k,t) = (p(—k,t). Thus, equation (|2.4|) is Hamiltonian and 
{(p(k,t),(p* (k,t) | k € N} are complex canonical variables. 

The Hamiltonian is conserved in time on smooth solutions, but it is 
an indefinite functional since it is cubically nonlinear. A positive-definite 
functional that is conserved on smooth solutions is the momentum 

V = Y^k(3*(k,t)(p(k) = — f \\d x \ 1/2 <p(x,t)] 2 dx. (2.8) 

The a priori estimate provided by the conservation of V (or the dissipation 
of V for weak solutions) is not sufficient to imply global existence, however, 
and we do not know of other conserved functionals of (|2.4|) for general T. 

2.2 Noncanonical variables 

It is convenient to introduce noncanonical dependent variables 
u:TxI ^R, m:ZxJ^C 

defined by 

u(x,t) = \d x \ 1/2 ip(x,t), u(k,t) = \k\ 1/2 fi(k,t). 
When rewritten in terms of n, equation (|2.4|) becomes 

oo 

ut(k,t)+ik S(—k,k — n,ri)u(k — n,t)u(n,t) = 0, (2.9) 

n=— oo 

where 

S(k,m,n) = ■ — for kmn ^ 0. (2.10) 

|/cmn| ' 

We define 5(0, m,n) = S(k,0,n) = S(k,m,0) = 0. The corresponding 
spatial form of (J2.9|) is then 

u t + a(u,u) x = 0, ( 2 -ll) 
where the bilinear form a is given by 

oo 

a(u,v)(k,t) = S(—k,k — n,n)u(k — n,t)v(n,t). (2.12) 
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The form a in (|2.11j) is related to b in H2.2|) by 

b{tp, iP) = \d x \ 1/2 a (\d x \ 1/2 ip, \d x \ 1/2 i?j . (2.13) 

It follows that equation ()2.11|) is invariant under the rescalings 

1 1— » rjt, t]x, u i ^ rf~^^ 2 u. 

As we will show in Section 14.21 under suitable assumptions on S the 
bilinear form a(u, v) may be regarded as a kind of 'nonlocal product' of u 
and v, so that ()2.11j) is a nonlocal generalization of the inviscid Burgers 
equation. Our main result is that if the kernel S(k, m, n) satisfies the bound 
in ()4.8|) below, where n = v — 3/2 > 0, and s > u, then the initial value 
problem for (|2.11|) - (|2.12|) has a unique solution for short times taking values 
in the Sobolev space H S (T). Equivalently, equation (|2.2|) - (|2.3|) has a short- 
time -fP-solution for s > v + 1/2. 

3 Examples 

In this section, we give a some examples of equations that belong to the class 
described in the previous section. The simplest, and standard, example is the 
inviscid Burgers equation, but this class of equations is surprising rich, and 
includes other interesting PDEs and integro-differential equations. These 
examples also illustrate various different possibilities for local and global 
existence. 

3.1 Burgers equation 

Suppose that 

T(k, m, n) = {kmnl 1 / 2 . (3.14) 

Then S(k, m,n) = 1 for kmn ^ 0, and the noncanonical equation (|2.11|l 
corresponding to (|3.14|) is the inviscid Burgers equation 

u t + (u 2 ) x = 0, (3.15) 

as can be seen from (|2.12|) by use of the convolution theorem. The kernel 
Q3. 14(1 is homogeneous of degree 3/2, corresponding to n = d in (|2.7|) . and 
this equation describes bulk waves, such as weakly nonlinear sound waves. 

Theorem ^ states local existence for (|3.15|) in H s with s > 3/2, which 
is the usual short-time existence result for first-order quasilinear hyperbolic 
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equations in one space dimension. Nontrivial periodic solutions form singu- 
larities in finite time, but the inviscid Burgers equation has global dissipative 
weak solutions that satisfy appropriate entropy conditions. These solutions 
are not, however, compatible with the Hamiltonian structure of (|3.15|) . 

3.2 Hunter-Saxton equation 

Consider the kernel 

T(k,m,n) = — Ikmnl 1 ^ 2 ( — H 1 ) for kmn ^ 0, (3.16) 

2 \ik im in ) 

with T(k, m,n) =0 when kmn = 0. Then 

S(k, m, n) = — ( — H 1 ) for kmn ^ 0. 

2 \ik im in J 

The corresponding noncanonical equation 1)2.11(1 is 

u t + (ud- l u) x -\{u 2 - (u 2 )) = 0, (3.17) 

where d~ l is the anti-differentiation operator on zero- mean, periodic func- 
tions, and the angular brackets denote the spatial average over a period. It 
follows that v = d~ 1 u satisfies the Hunter-Saxton equation |14j . 

vt + vv x - l -d- x [yl - (v%)) = 0. (3.18) 

The kernel (|3.16|) has degree v = 1/2, corresponding to n = d — 2 in 
l|2.7jl . This is consistent with the derivation of (|3.18|) in Jl] for orientation 
waves in a massive director field. The waves propagate in three-dimensional 
space, so d = 3, and the mass parameter of the medium is set by the moment 
of inertia per unit volume of the director field, which has dimension M/L, 
so n = 1. 

Theorem |21 (with A = 1) implies the short-time existence of IP-solutions 
of (l3T7l for s > 1/2, and hence iP-solutions of (ETTSl for s > 3/2. This 
agrees with the short-time existence result of Yin [S^. Nontrivial periodic 
solutions of the Hunter-Saxton equation form singularities in finite time, 
but the equation has global weak solutions, including both dissipative and 
energy-conserving weak solutions (|S], P4 ). The conservative weak 

solutions are compatible with the Hamiltonian structure of the equation. 
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It is interesting to note that (|3.18|) is completely integrable (|I], |16j). 
and also arises as the high-frequency limit of the integrable Camassa-Holm 
equation [7j, 

v t + 2kv x + 3vv x = (v t + w x ) xx - - (v 2 x ) x . 
3.3 Surface wave equations 

As discussed in the introduction, weakly nonlinear asymptotics for 'genuine' 
hyperbolic surface waves typically leads to nonlocal scale-invariant equa- 
tions. These equations inherit a Hamiltonian structure from the Hamilto- 
nian structure of the primitive physical equations (such as hyperelasticity 
or magnetohydrodynamics). 

Perhaps the simplest kernel arising for surface waves is given by 

T(k,m,n)= | ,' l | mn | | , , (3.19) 
\k\ + \m\ + \n\ 

leading to the spectral evolution equation 

E2 1 /l 1 1 /l T\i 1 1 T\ I 
,,, : |, j— r -Mk - n,t)(p(n,t) = 0. (3.20) 
fc + \k — n + \n\ 

n=—oc 111 iii 

This equation was proposed by Hamilton et al ( jlUj. |llj ) as a model equa- 
tion for nonlinear Rayleigh waves, and it describes the propagation of surface 
waves on a tangential discontinuity in incompressible magnetohydrodynam- 
ics p. 

The kernel in (|3.19|) is homogeneous of degree v = 2. This degree is 
consistent with the result of dimensional analysis in (|2.7|) . since the dimen- 
sion d of the wavenumber space of surface waves that propagate along the 
boundary is one less that the spatial dimension n of the medium (whose 
density sets the mass parameter of the wave motion). 

The spatial form of (|3.2()|) may be written explicitly as (|2j, jllj ) 



ipt + \n[^ 2 ] xx + ^ xx = 0, (3.21) 



2 

where H is the Hilbert transform and ip = H[yj]. This is equation (jl.lj) . 
Equation (|3.21j) corresponds to ()2.2j) with 

Hf,f) = \ {^ 2 ) xx -n[^xx] 

3 f 

= J ^ip x <f x dx. 
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The appearance of second-order spatial derivatives in the expression for 
b is misleading since some cancelation occurs. It follows from (j2.13|) and 
Proposition [T] (with \i = 1/2) that b : H s x H s — > H^ 1 involves the 'loss' 
of only one derivative when s > 3/2. 

Theorem ^ implies that the noncanonical equation associated with (|3.21|) 
has short-time iP-solutions when s > 2, so Q3.21JI has .fP-solutions when 



Numerical computations (pQ, [57]) show that solutions of ()3.21j) 

form singularities in which the derivative <p x blows up, but (p appears to 
remain continuous. The global existence of appropriate weak solutions is an 
open question. As far as we have been able to determine, equation (|3.21|) is 
not completely integrable. 

Similar equations describe the evolution of weakly nonlinear Rayleigh 
waves on an elastic half-space. In that case, <p(x, t) corresponds to the ver- 
tical displacement of the boundary of the half-space. For weakly nonlinear 
Rayleigh waves on an isotropic elastic half-space, cp(k, t) satisfies an equa- 
tion of the form with the kernel (see Qlj, [H], EB, HH , US, |S2], for 
example) 



Here, < r = ct/ci < 1 is the ratio of the elastic solid's transverse and 
longitudinal wave speeds, q and q respectively, and a, 0, 7 are nonlinear 
elastic constants. The corresponding equation for (p(x, t) has short-time 
-£P-solutions for s > 5/2. 

More complicated, but qualitatively similar, kernels arise for nonlinear 
Rayleigh waves on a non-isotropic elastic half-space ( ^2] > H2j ) • Kernels with 
a different degree of homogeneity (y = 5/2) arise in the weakly nonlinear 
description of elastic edge waves ([213]) |2S])- 

3.4 Compacton equations 

The kernel 



s > 5/2. 



T(k, n, m) 




(3.22) 



T(k, m,n) = i \kmn 



2 kmn 



(3.23) 
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leads to a noncanonical equation Q2.11JI that is a local PDE, 

n * + = °- 

Differentiating this equation with respect to x and setting u x = v, we get 

* + (A™ = °- ( 3 - 24 ) 

This equation is the high-frequency limit of the K(2, 2)-equation introduced 
by Rosenau and Hyman |28j : 

vt + {v 2 ) x + (v 2 ) xxx = 0. (3.25) 

The K(2, 2)-equation has 'compacton' solutions (that is, compactly sup- 
ported traveling waves) given by 

4c ( x ~ 
v(x, t) = — cos I — - — J if \x — ct\ < 2tt, 

and v(x, t) = if \x — ct\ > 2ir. Some numerical solutions of IJH.25J) are given 
in 0, |23| . 

Although (|3.24|) does not have compacton solutions, it has travelling 
wave solutions that describe the local non-smooth behavior of the compacton 
solutions of (|3.25|) near the 'edges' of the compacton; for example: 

v(x, t) = ^ (x — ct) 2 if x < ct, 

and v(x, t) = if x > ct. The relationship between (|3,24l) and (|3.25|) is thus 
similar to the one between the Hunter-Saxton and Camassa-Holm equations, 
including the local description of the non-smooth behavior of the 'peakon' 
solutions of the Camassa-Holm equation by solutions of the Hunter-Saxton 
equation. As far as we know, however, neither Q3.24JI nor ()3.25j) is completely 
integrable. 

An expansion of the spatial derivatives in ()3.24|) , 

v t + 2vv xxx + 6v x v xx = 0, 

leads to a backward diffusion term when v x > 0, which suggests that the 
equation may not be well-posed. The short-time existence proof given here 
does not apply to (|3.24|) since its kernel does not satisfy ()4.9|) below, and 
this condition seems to be necessary to obtain L 2 -Sobolev estimates. It is 
conceivable that other estimates could lead to short-time existence for (|3.24|) 
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or (|3,25|) . but no such results appear to be known (although it is possible 
that a Cauchy-Kowalevski type argument could be used for analytic initial 
data 0). 

More generally, the kernel 

T(k,m,n) = |/cmra| 1//2 {ikmnf 

of degree v = 3p + 3/2, has the corresponding noncanonical equation 

Ut + d p x +l {d p x uf = o, 

which may be regarded as a higher-order generalization of the inviscid Burg- 
ers equation, to which it reduces whenp = 0. Our short-time existence result 
does not apply to this equation when p > 0. 



4 Local existence of smooth solutions 

We consider an initial value problem for the noncanonical form of the equa- 
tion described in Section |2J 

u t + a(u,u) x = 0, (4.1) 
u(x,0) = f{x). 

Here, u : T x I — > R, with I an open time interval containing the origin, 
/ : T — > R is given initial data, and the bilinear form a is defined by 

oo 

a(u,v)(k,t) = S(—k,k — n,n)u(k — n,t)v(n,t). (4.2) 

n=— oo 

We assume that the kernel S : 1? — > C satisfies 

S(k,m,n) = if kmn = 0, (4-3) 

S(k,m,n) = S* \—k, —m, —n), (4-4) 

S(k,m,n) = S(k,n,m), (4.5) 

S(k,m,n) = S(m,n,k), (4-6) 

S(rjk,r]m,r]n) = if S(k,m,n) for 77 > 0, (4-7) 

These properties follow from those of T in (|2.5j) and the definition of S in 
(PHUjt . with 

3 

\i = v . 
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We will focus our attention on kernels S with nonnegative degree fi > 0. We 
consider kernels of negative degree in Section [4,41 

In order to develop a local existence theory, we assume that S satisfies 
the following condition for some constant C and all k,m,n € Z: 

\S(k,m,n)\ < Cmxn{\k\^,\mf,\n\^}) (4.8) 

equivalently, T satisfies 

\T(k,m,n)\ < C \kmn\ 1/2 min [\k\ u ~ 3 / 2 , |m|^ 3 / 2 , |n|^ 3 / 2 } . (4.9) 

For example, the kernel ()3.14j) for the inviscid Burgers equation satisfies 
with v = 3/2 and C = 1. The surface- wave kernel ({SUSt satisfies (|43|) 
with ^ = 2 and (7=1, since 

2|fcmn| 1 / 2 . f.-n/o , .i/o , ,1/2 1 

,,, 1 1 , , , < mm <^ (/el 1 / 2 , (ml 1 / 2 , \n\ 1 ^ \ , 
\k\ + jm| + \n\ (. j 

and the Rayleigh-wave kernel (|3.22|) satisfies ()4.9|) with v = 2 and 

C= i/3H +M + 



2 

By contrast, the compacton-type kernel (|3,23|) does not satisfy ()4.9j) . 
4.1 Notation 

We denote by H S (T), or .fP for short, the Hilbert space of periodic functions 
with zero mean and square integrable derivatives of the order s, 

H s (T) = h:T^R\ jr \k\ 2s \f(k)\ 2 <oo, /(0) = 1 . 

I. k=— 00 J 

As inner product and norm on H s , we use 



(f,g) s = J2 \k\ 2s f(k)g(-k) 

k=—oo 

11/11, = ( E \k\ 2s f(k) 2 ] 



\k=— 00 
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In particular, 

For / : Z — > C, we also define 

|.f 

so that 

4.2 The bilinear form 



E |/(*) 



1/2 



1/p 



\k=— oo 



l/ll 



1*17 



f 2 



(4.10) 



First, we obtain some properties of the bilinear form a given by (|4.2[) . In 
the simplest case when m, n) = 1 for fcmn 7^ 0, corresponding to /i = 0, 
the form a is a multiplication operator on H s for s > 1/2: 



a(u, v) = uv 



1 

27 



uv dx. 



(4.11) 



As the following proposition shows, the form a has similar properties to a 
multiplication operator for any kernel S that satisfies (|4.8|) . although we 
need some additional smoothness to ensure that .£P is an algebra under the 
product defined by a. 

Proposition 1 Let a(u, v) be defined by \4-°>$) , where S : Z 3 — > C satisfies 
\4.'J!j) - \4.tfj) . Suppose that s > fx + 1/2. 

(a,) TTten a : iP x iP — > iP is a bounded symmetric bilinear form, 
(b) For all u,v G iP+\ 



a(u,v) x = a(u x ,v) + a(u,v x ). 
(c) For all u,v,w € iP, 

1 



2^ 



— / ua(v,w)dx = — / va(w,u)dx 



2tt 



and i/iere is a constant M s _^ > suc/i i/iai 



— / ua(v,w)dx 



< CM s -,j,\\u\\b\\v\\o\\w\\o. 



(4.12) 



(4.13) 



(4.14) 
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Proof. Using (|4.2j) . we get 



E 



|A;| S S(—k, k — n, n)u{k — n)v(n) 



(4.15) 



For s > and m, n € M, we have 

|m + n| s < B s (\m\ s + |n| 

where 



1 ifO<s<l, 
2 s - 1 if s > 1. 



Using this inequality, with m = k — n, and (|4,8|) in (|4.15j> . we get 

||o(Ti,«)||g < ^ X] ^ {\k - n\ s + [n| s ) \S(-k,k- n,n)u(k - n)v{n)\ 

k n 

< B 2 S C 2 5^ (I* - H s |n|^ + | A; - n|^|n| s ) \u(k - n)v{n)\ 



< B 2 s C 2 \\(\k\ s \u\) * {\kf\v\) + * (\k\ s \v\) 

Here, we denote the convolution of f,g : Z — > C by 

oo 

/*?(*) = f( k - n )a(n)- 



It follows that 

\\a(u,v)\\ s < B S C [\\(\k\ s \u\) * (\k\"\v\) 



+ 



(ifenfii)*^!-!^)^}. 



Using Young's inequality, 



f*9 



< 



t 2 



c 1 



(4.16) 



we get that 

\\a(u,v)\\ s < B s c{\\u\ 



e 1 



+ 



\k\"\u\ J|«||,}. (4.17) 
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If s > 1/2, we have for all / G H s that 



/ 



c 1 



<M a \\f\\ 



where 



' oo % 
2 n 2s 



1/2 



(4.18) 
(4.19) 



71=1 



Hence, if s > // + 1/2, inequalities (|4.17|) and Q4.18JI imply that 

||a(«,u)|| s < 2B s CM s - fl \\u\\ s \\v\\ s . 

Thus, a is bounded on H s . The symmetry and bilinearity of a is obvious. 
To prove (|4.12|) . we observe that 

a(u,v) x (k) = ika{u,v)(k) 

= ik S( — k, k — n, n)u{k — n)v(n) 
n 

= i(k — n)S(—k, k — n, n)u{k — n)v(n) 

n 

+ inS(—k, k — n, n)u(k — n)v(n) 

n 

= a(u x ,v)(k) + a(u,v x )(k). 
Equation Q4.13JI follows from the fact that 

— / ua(v,w)dx = u(—k)a(v,w)(k) 
2vr J T ^ 

= S(—k, k — n, n)u{—k)v{k — n)w(n) 



k,n 



is symmetric under cyclic permutations of u, v,w. Moreover, using (|4.8jl . the 
Cauchy- Schwartz inequality, Young's inequality, and (|4.18|) in this equation, 
we get for s > fj, + 1/2 that 



1 

2^ 



ua(v, w) dx 



< C^2\k\^\u{-k)v(k-n)w{n)\ 







fe,n 


< 


chk\i>\ 


< 


c 


\k\»u 


< 


| 



[10 \ 



p 

\v\\p \\w\\p 



p 



□ 
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The form a satisfies other Sobolev inequalities. In Proposition ^ we 
prove only the ones that we need to use below. 



4.3 Local existence 

We define 

K s = CC s M s ^x, (4.20) 

where C is the constant in (|4,8|) . C s is the constant in ()5.2j) . and M s is given 
by (l4~m 

Theorem 1 Suppose that S : Z 3 — > C satisfies j4-^ ^ ^4~ 7 ^ an d \4-$[) with 
fi > 0, the form a is given by \4-^ , and s > /i + 3/2. Then, for any 
f £ H S (T), the initial-value problem j4-l\ ) has a unique local solution 

uecfl; H S (T)) n C 1 (7; H S -\T)) (4.21) 



defined on the time interval I = (— £*,£*), where 



Proof. Let P N : H S (T) — > H S (T) denote the orthogonal projection 

/ oo \ N 



p n [ e f( k y kx ) = E f( k y kx - ( 4 - 23 ) 

Vfc=-oo / k=-N 

We introduce the Galerkin approximations 

N 

u N (x,t) = E u N (k,t)e ikx 

k=-N 

that satisfy the system of ODEs 

u? + P N a(u N ,u N ) x = 0, (4.24) 
u N (x,0) = P N f(x). 

The main estimate for the local existence proof is given in the following 
Proposition. 



18 



Proposition 2 Suppose u N (x,t) satisfies 



d 

-r \\U 

dt 



N 



< K< \\u 



TV 1 1 2 



and s > fj, + 3/2. Then 

(4.25) 



where K s is defined in 



Proof. Dropping the superscript N to simplify the notation, we compute 
from flOjj) , (jOHjl . and (JO) that 

|^|fc| 2 ^(M)£(-M) 

+2i ^ fc|Jfe| 2a 5 , (-Jfe, /c - n, n)u(-k, t)u(k - n, t)u(n, t) = 0. (4.26) 
fc,n 

Using the cyclic symmetry of S(k,m,n), we can rewrite this equation as 
d 



2a 



n)|fe - n| 2s - n|n| 2s } 



fc.n 



S(—k, k — n, n)u(—k, t)u{k — n, t)u(n, t) = 0. 
Using the inequality l|5.2jl . we find that 
d 



dt 



E^i 28 \u(k,t)[ 



2C 



< — -^{|A;| s |fc-n|>| + jA;-n| s |n| s |A:| + \n\ s \k\ s \k - n\] 



\S(—k,k — n,n)| \u(— k,t)u(k — n,t)u(n, t)\ . (4.27) 

Using (|4,8|) in the right hand side of (|4.27l) . and applying the appropriate 
bound on each term, we get 



£>| 2s RM)f 



2CC, 



< — — ^{\k\ s \k - n\ s \n\» +1 + \k- n\ s \n\ s \k\^ +1 



kji 



+ |n| s |A:| s |A; - n|' t+1 } \u(-k,t)u(k - n,t)u(n,t)\ 
< 2CC s ^\k\ s \k - n\ s \n\» +1 \u(-k,t)u(k - n,t)u{n,t)\ 



fc.n 
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Using the Cauchy-Schwartz inequality and Young's inequality (|4,16|) on the 
right hand side of this inequality, we get 



d 




2 


dt 


\k\ s u 


p 



< 2CC S 


\k\ s u 


i 2 






2 


< 2CC S 


\k\ s u 


I 2 



It then follows from (j4~lT)|) and (|4~T%|) that if s > fj, + 3/2, then 



d_ 

It 



\u\ 



< 2CC S M S 



-1 u 



Simplifying this equation and using (|4.20j) we get (|4.25|) . □ 
Using Gronwall's inequality, we deduce from (|4.25|) that 

h N Ut)< 



1-K S 



\t\' 



so the Galerkin approximations exist in \t\ < i*, where is given by (|4.22|) . 

It follows from Proposition ^ an d standard arguments (for example, |3(J| . 
§16.1, or ^Hl) that we can extract a subsequence (u ) of Galerkin approxi- 
mations such that 



,iV 



ueC w [J;H 



(4.28) 



on any closed interval J C I, where u is a distributional solution of (|4,1|) . 
(Here, C w denotes the space of weakly continuous functions.) Proposition |^1 
and Gronwall's inequality imply that t i— > ||n(t)|| s is continuous. It follows 
that u is strongly if s -continuous, and thus satisfies (|4.21|) since a is strongly 
continuous on H s . 

To prove uniqueness, we suppose that u, v are two solutions of ()4.1|) 
with the regularity stated in (|4.21j) . Then, using integration by parts and 
the properties of a given in Proposition ^ we compute that 



d ,, ,,2 

— \\U — V\\n 

dt 11 110 



— / (u — v) 2 dx 

2tt dtJ T y ' 



- [ (u-v)[a(u,u) x -a(v,v) x ] dx 



1 

7T Jt 



(u — u) a(u + v,u — v) x dx 
(u — v) x a(u + v, u — v ) dx 
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— / (u + v) a (u — v, u x — v x ) dx 

— j (u + v) a (u — v,u — v) x dx 
1 

~2tt 



T 



r 



{u + v) x a(u — v,u — v) dx. 



Since s > /x + 3/2, it follows from this equation and (|4.14j) that 



d 



— \\u — v\ 



< CM s _ M _i||n + v\\ s \\u - v\\l 



dt 

Hence, u = v by GronwalFs inequality, and the solution is unique. □ 
4.4 Negative degree kernels 

In this section, we consider kernels S of negative degree /i = —A < that 
satisfy the following condition for some constant C and all kmn ^ 0: 



\S(k,m,n)\ < ^ nux (4.29) 



C_ 

min |m| A , 

This corresponds to the assumption that T satisfies the condition 



C Ikmnl 1 / 2 

\T(k,m,n)\ < r- — rrr^ ; — tttt; ; — tttt; — r- (4.30) 

n ~ min{\k\ 3 / 2 - v ,\m\ 3 / 2 -' / ,\n\ 3 / 2 - u } v ' 

where v = 3/2 — A < 3/2. In particular, Q4.3UJ1 is satisfied by the Hunter- 
Saxton kernel (j3~Tp|) with v = 1/2 and C = 3/2. 
Let 

K s ,\ = CC Sj \M s+ \-x, 

where C is the constant in (|4.29|1 . C s ^\ is the constant in ()5.1j) . and M s is 
given by (|4TP3j) . 

Theorem 2 Suppose that S :1? -> C satisfies £pfy-M7fy and with 
A > 0, the form a is given by \4-°4) j an d 

/ 3 x 1 

s > max <^ - - A, - 

Then, for any f £ H S (T), the initial-value problem j4.1\) has a unique local 
solution 

ueC(l; H S (T)) n C 1 (i; iT s-1 (T) N 
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defined on the time interval I = (— where 

1 



K. 



Proof. The proof is essentially the same as before. For s > 1/2, we get that 
a : H s x H s — » i/ s is a continuous bilinear form. From (|4.26l) and (|4.29|) . 
the Galerkin approximations satisfy 



d 



-Y^^ \u(k,t) 



< 



J2\k\k\' 



n)\k — n\ 2s — n\n 



k,n 



\S(—k, k — n, n)\ \u(— k, t)u(k — n, t)u(n, t)\ 



^ 2C ^ |fc|fc| 2 * -(k-n)\k- n| 2s - n\n\ 2s \ 

- ~T~ ■ m ia \ U ia i \\x \u{-k,t)u{k - n,t)u{n,t)\ 

S mm^ A; , \k — n , \n A \ 

k,n 

Using (|5.1|) in this equation, and summing over cyclic permutations, we find 
that 



J2\k\ 2s \u(k,t)f 



< 2CC s>x J2\ k \ S \ k 



n\ n 



l-A 



k,n 



< 2CC, 



s.A 



u(—k, t)u(k — n, t)u(n, t) | 
\k\ s u 2 



Ikl^u 



c 1 



p 



Hence, using (|4.18|) . with s > 3/2 — A, and simplifying the result, we get 

d „ 



< CC S , A M S+A _ 1 ||u| 
and the local existence result follows. □ 



5 An Inequality 

For completeness, we prove an inequality used in our local existence proof. 
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Proposition 3 For s > and A > 0, there exists a constant C Sj a such that 
for all nonzero k, m, n £ R to^/j /c + m + n = 0, 



|/c|/j| 2s + m|m| 2s + n|n 



2.s 



min |m| A , 

< C^ll^H^nl 1 ^ + Inn^nml 1 ^ + HX 5 ^ 1 ^} . (5.1) 

Proof. We can assume without loss of generality that < m < n and 
A; = — [m + n) < 0. Writing x = m/n, we have 



l^l/cl 2,5 + m\m\ 2s + n|n.| 2s | 



|m| A {\k\ s \m\ s \n\ l A + |n| s | A;| s |m| 1 A + |ttz.| s |n | s | | 1 A } 
where / S)A : (0, 1] -> (0, oo) is given by 

(x + l) 2s+1 - x 2s+1 - 1 



fs,\{x), 



fs,x(x) 



(x + l) s (x s+x + x) + (x + l) 1 - A x s+A ' 

This function is continuous on (0, 1] and has a finite limit as x — > + , with 

[ if < s + A < 1, 

/.,A (0 + ) = < (2s + l)/3 if s + A = 1, 
[ 2s + 1 if s + A > 1. 

Hence / S) ^ is bounded on (0, 1] by C s> \, say, which proves ()5.1|) . □ 

In particular, if A = 0, we write f s = f s> o and C s = C s ,o, so that 

|/c|/c| 2s + m\m\ 2s + n|n| 2s | 

< C s {|fc| s |m| s |n| + \n\ s \k\ s \m\ + |m| s |n| s |A;|} (5.2) 

for all k, m, n G R with + m + n = 0. We have fi(x) = 1, fe{x) = 7, and 
f s (l) = 2 s — 1. Numerical plots show that / s is monotone decreasing on 
(0, 1] for 1 < s < 3 and monotone increasing for s > 3, in which case 

/ 2s + 1 if 1< s < 3, 
Cs ~\ 2 s -1 ifs>3. ^ 



23 



References 



[1] G. Ali, and J. K. Hunter, Nonlinear surface waves on a tangential dis- 
continuity in magnetohydrodynamics, Quart. Appl. Math., 61 (2003), 
451-474. 

[2] G. Ali, J. K. Hunter, and D. F. Parker, Hamiltonian equations for 
scale-invariant waves, Stud. Appl. Math., 108 (2002), 305-321. 

[3] M. Artola, and A. J. Majda, Nonlinear development of instabilities in 
supersonic vortex sheets. I. The basic kink modes, Physica D, 28 (1987), 
253-281. 

[4] R. Beals, D. H. Sattinger, and J. Szmigielski, Inverse scattering solu- 
tions of the Hunter-Saxton equation, Appl. Anal. 78 (2001), 255-269. 

[5] S. Benzoni-Gavage, F. Rousset, D. Serre, and K. Zumbrun, Generic 
types and transitions in hyperbolic initial-boundary-value problems, 
Proc. Roy. Soc. Edinburgh Sect. A 132 (2002), 1073-1104. 

[6] A. Bressan and A. Constantin, Global solutions of the Hunter-Saxton 
equation, preprint. 

[7] R. Camassa and D. Holm, An integrable shallow water equation with 
peaked solitons, Phys. Rev. Lett. 71 (1993), 1661-1664. 

[8] O. Costin and S. Tanveer, Existence and uniqueness for a class of higher- 
order partial differential equations in the complex plane, Comm. Pure. 
Appl. Math. 53 (2000), 1092-1117. 

[9] J. de Frutos, M. A. Lopez-Marcos, and J. M. Sanz-Serna, A finite dif- 
ference sceme for the K(2, 2) compacton equation, J. Comput. Physics 
120 (1995), 248-252. 

[10] M. F. Hamilton, Yu. A. Il'insky, and E. A. Zabolotskaya, Local and 
nonlocal nonlinearity in Rayleigh waves, J. Acoust. Soc. Am., 97 (1995), 
882-890. 

[11] M. F. Hamilton, Yu. A. Il'insky, and E. A. Zabolotskaya, Evolution 
equations for nonlinear Rayleigh waves, J. Acoust. Soc. Am., 97 (1995), 
891-897. 

[12] M. F. Hamilton, Yu. A. Il'insky, and E. A. Zabolotskaya, Nonlinear 
surface acoustic waves in crystals, J. Acoust. Soc. Am., 105 2 (1999), 
639-651. 



24 



[13] J. K. Hunter, Nonlinear surface waves, Contemporary Mathematics, 
100 (1989), 185-202. 

[14] J. K. Hunter and R. Saxton, Dynamics of director fields, SIAM J. Appl. 
Math. 51 (1991), 1498-1521. 

[15] J. K. Hunter and Y. Zheng, On a nonlinear hyperbolic variational equa- 
tion. I. Global existence of weak solutions, Arch. Rational Mech. Anal. 
129 (1995), 305-353. 

[16] J. K. Hunter and Y. Zheng, On a completely integrable nonlinear hy- 
perbolic variational equation, Phys. D 79 (1994), 361-386. 

[17] N. Kalyanasundaram, Nonlinear surface acoustic waves on an isotropic 
solid, Int. J. Eng. Sci., 19 (1981), 279-286. 

[18] N. Kalyanasundaram, R. Ravindran and P. Prasad, Coupled amplitude 
theory of nonlinear surface acoustic waves, J. Acoust. Soc. Am., 72 
(1982), 488-493. 

[19] T. Kato, Nonlinear equations of evolution in Banach spaces, Proc. 
Symp. Pure Math., 45 (1986), 9-23. 

[20] A. P. Kiselev, Nondispersive edge waves, Wave Motion, 29 (1999), 111— 
117. 

[21] R. W. Lardner, Nonlinear surface waves on an elastic solid, Int. J. 
Engng. Sci., 21 (1983), 1331-1342. 

[22] R. W. Lardner, Nonlinear surface acoustic waves on an elastic solid of 
general anisotropy, J. Elast., 16 (1986), 63-73. 

[23] D. Levy, C-W Shu, and J. Yan, Local discontinuous Galerkin methods 
for nonlinear dispersive equations, J. Comput. Phys. 196 (2004), 751- 
772. 

[24] A. P. Mayer, Surface acoustic waves in nonlinear elastic media, Phys. 
Rep. 256 (1995), 237-366. 

[25] D. F. Parker, Waveform evolution for nonlinear surface acoustic waves, 
Int. J. Engng. Sci., 26 (1988), 59-75. 

[26] D. F. Parker, Elastic edge waves, J. Mech. Phys. Solids, 40 (1992), 
1583-1593. 



25 



[27] D. F. Parker, and F. M. Talbot, Analysis and computation for nonlinear 
elastic surface waves of permanent form, J. Elast., 15 (1985), 389-426. 

[28] P. Rosenau and M. Hyman, Compactons: solitons with finite wave- 
length, Phys. Rev. Lett. 70 (1993), 564-567. 

[29] M. Sable- Tougeron, Existence pour un probleme de Pelastodynamique 
Neumann non lineaire en dimension 2, Arch. Rational Mech. Anal. 101 
(1988), 261-292. 

[30] M. Taylor, Partial differential equations III Nonlinear equations, 
Springer- Verlag, New York, 1997. 

[31] Z. Yin, On the structure of solutions to the periodic Hunter-Saxton 
equation, SI AM J. Math. Anal. 36 (2004), 272-283. 

[32] E. A. Zabolotskaya, Nonlinear propagation of plane and circular 
Rayleigh waves in isotropic solids, J. Acoust. Soc. Am. 91 (1992), 2569- 
2575. 

[33] V. E. Zakharov, V. S. L'vov, and G. Falkovich, Kolmogorov Spectra of 
Turbulence I, Springer- Verlag, Berlin, 1992. 

[34] P. Zhang and Y. Zheng, On the existence and uniqueness of solutions 
to an asymptotic equation of a variational wave equation, Acta Math. 
Sinica 15 (1999), 115-130. 



26 



